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We examine the prospects of probing nonstandard interactions (NSI) of neutrinos in the e — T sector with 
upcoming long-baseline — > Vg oscillation experiments. First conjectured decades ago, neutrino NSI remain 
of great interest, especially in light of the recent solar neutrino measurements by SNO, Super-Kamiokande, 
and Borexino. We observe that the recent discovery of large 613 implies that long-baseline experiments have 
considerable NSI sensitivity, thanks to the interference of the standard and new physics conversion amplitudes. 
In particular, in some parts of NSI parameter space, the upcoming NOvA experiment will be sensitive enough 
to see ~ 3(7 deviations from the SM-only hypothesis. On the flip side, NSI introduce important ambiguities 
in interpreting NOvA results as measurements of CP-violation, the mass hierarchy and the octant of 023- In 
particular, observed CP violation could be due to a phase coming from NSI, rather than the vacuum Hamiltonian. 
The proposed LBNE experiment, with its longer --^ 1300 km baseline, may break many of these interpretative 
degeneracies. 
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I. INTRODUCTION 

"The effect of coherent forward scattering must 
be taken into account when considering the oscil- 
lations of neutrinos traveling through matter. In 
particular [. ..] oscillations can occur in matter 
if the neutral current has an off-diagonal piece 
connecting different neutrino types. Applications 
discussed are solar neutrinos and a proposed 
experiment involving transmission of neutrinos 
through 1000 km of rock." 

Though the above quote could easily have been written 
this year, or even applied to the present paper, it was writ- 
ten presciently in 1978 by Lincoln Wolfenstein in his semi- 
nal paper on the effects of matter on neutrino oscilations ||T|. 
Although originally proposed as an alternative to mass in- 
duced oscillations OJBl, beyond-the-Standard-Model (BSM) 
neutrino-quark interactions remain a phenomenological possi- 
bility (e.g., IIS HTTll ) that can produce potentially observable ef- 
fects in oscillation experiments. Three decades after the above 
quote was written, we have finally reached the era of 1000 km 
experiments, with several years of data collected at MINOS, 
NOvA launching next year, and LBNE on the drawing board. 

Our goal in this paper is to gauge the sensitivity of these 
experiments to NSI, in light of what has become known 
about neutrino oscillations over the last decade. We delib- 
erately choose to avoid a full analysis that scans over many 
couplings with different flavor combinations and consider 
a simplified framework with only one effective flavor off- 
diagonal piece connecting electron- and tau-type neutrinos, 
^ D -2V2 e/rGf (fy^ifVeY'^VT:) +h.c., where / = u,d,e. 
We will see that this framework nonetheless reveals a rich 
spectrum of physical possibilities. Importantly, Eer has its own 



'Electronic address: friedland@lanl.goyl 
^Electronic address: ianshoe@lanl.gov 



CP-violating phase and can lead to ambiguity in interpreting 
the searches of Cf -violation and the mass hierarchy. 

As a first illustration, let us examine the effect this one pa- 
rameter can have on the solar electron neutrino survival prob- 
ability, P{Ve — ^ Ve). The standard large mixing angle (LMA) 
MSW solution makes a definite prediction for how this prob- 
ability varies as a function of the neutrino energy, Ey. This 
prediction is shown as a dashed line in Fig. [T] with the gray 
band around it coming from the uncertainty on the standard 
oscillation parameters. Both are taken from |18 |. Also taken 
from 1 18 1 is the allowed region of this probability inferred 
from all three stages of SNO data, as labeled in Fig. [T] At 
low energies, we also include the survival probabiUties of pp 
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FIG. 1 : Recent SNO solar neutrino data |T8| on P{Ve — > Ve) (blue line 
with 1 (7 band). The LMA MSW solution (dashed black curve with 
gray 1 CT band) appears divergent around a few MeV, whereas for 
NSI with £er = 0.4 (thick magenta), the electron neutrino probability 
appears to fit the data better The data points come from the recent 
Borexino paper 1191 . 
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and ''Be neutrinos and the pep neutrinos as recently reported 
by Borexino |19]. While the standard MSW solution pro- 
vides an acceptable fit to all this data, the addition of the gg^ 
NSI parameter noticeably improves the fit, as shown with the 
thick magenta curve. This finding is consistent with earlier 
work |l9l[20l, although notice that here we vary only the off- 
diagonal e — T coupling, consistent with our framework. Note 
that although we have chosen Eez = 0.4 in Fig.[T](see II B for 
the definition of Eer), there is a family of continuous defor- 
mations that approach the MSW LMA solution as Eex goes to 
zero. The NSI analysis in |20|, which preceded the three-stage 
SNO data release 1 18 1, found a best fit of £ex — 0.2. We hence 
focus on \£i,t\ of the order 0.2 — 0.4 as fiducial values. 

Importantly, this choice is consistent with the wealth of 
other experimental information. For example, atmospheric 
neutrino data from Super-Kamiokande allow values of \£ex\ ^ 
0.5 lfT0l - [T2l [TtII . and become even less restrictive if tunings 
between different NSI parameters fTTl are allowed. A good 
fit to the data with nonzero Eex is achieved by shifting the vac- 
uum oscillation parameters 023 and A32 (see the definitions 
below) as shown in Fig 2 of fTTI. Since the direction of the 
required shift roughly coincides with the direction of reduced 
sensitivity at MINOS 1 12 1, even despite recent MINOS data, 
large values of NSI remain allowed 1 17 |. 

Similarly, the bounds from non-oscillation experiments do 
not exclude values of Eex hinted at by solar data. For exam- 
ple, quark-NSI couplings at hadron colliders manifest them- 
selves as anomalous monojet, monophoton II2TII . and multilep- 
ton signals 1211 l22l . Presently, the LHC and the Tevatron are 
starting to encroach on interesting regions of parameter space 
for heavy or resonantly produced mediators, but provide es- 
sentially no constraint when the mediator mass is < 10 GeV. 
(Similar constraints on NSI involving electrons were derived 
from LEP data using monophoton events ll23l .) It is worth em- 
phasizing that this "light mediator" regime could easily first 
appear in solar neutrino or long-baseline data. A similar pos- 
sibility of very light (<C GeV) mediators in the neutrino sector 
has been recently invoked in order to explain dark matter di- 
rect detection anomalies ll24H26ll . the short-baseline anoma- 
lies 1271 , and the appai-ent CPT violation at MINOS ||28| 
(which has since been resolved). 

Other constraints on NSI in the contact operator limit, from 
fixed target experiments and rare decays, are compiled in |8|. 
For example, for the Eex coupling in question, the CHARM 
beam-dump experiment |29| constrains Eg^ < 0.5, while con- 
straints on hadronic decays of the tau restrict ^ 1-6 EJ. 
A precise determination of how such experimental constraints 
are altered when the mediator of neutrino NSI is allowed to 
be light is yet to be carried out. Moreover, light mediators 
of NSI can contribute to quark-quark and neutrino-neutrino 
scattering Stringent bounds on light vector mediators between 
quarks exist from low energy neutron-nucleus scattering ex- 
periments 1 30 1, while supernovae may provide a constraint on 
neutrino-neutrino scattering. The analysis becomes model- 
dependent and is beyond the scope of this work. 

The further testing of the anomalous matter effect due to 
NSI should occur with the next generation of solar neutrino 
experiments, notably SNO-n, and with the long-baseline ex- 



periments, notably NOvA. Both are currently under con- 
struction and are expected to start operations next year Af- 
ter NOvA, the LBNE experiment, with a longer 1300 km 
baseline is expected to come online. The NSI sensitivity of 
NOvA and LBNE is the subject of the present paper. As we 
show below, these experiments may probe NSI couplings at 
the levels hinted at by the solar data. This is made possi- 
ble by the recently measured large value of 613, such that the 
dominant NSI contribution comes at linear order, from in- 
terference between standard and non-standard physics. These 
searches will be especially of great interest should SNOh- con- 
firm the distortion of the ^B solar neutrino spectrum. 

The outline of this paper is as follows. In Sec. |II A| we 
specify our conventions for the description of standard oscil- 
lation physics. In Sec. II B we introduce the formalism for 



describing NSI and derive a useful analytic expression for the 
Ve conversion probability, highlighting that the domi- 
nant effect from NSI comes from interference between stan- 



dard and non-standard physics. In Sec. Ill we show that MI- 
NOS is at the edge of sensitivity to a portion of the parameter 
space favored by solar data, though with significant degen- 



eracies of interpretations existing. In Sec. IV we turn to the 



NOvA experiment and show that a combined analysis of the 
neutrino and antineutrino conversion probabilities may reveal 
a measurable effect of the NSI. We also discuss how the pos- 
sible presence of NSI degrades NOvA's ability to determine 
the hierarchy and the vacuum phase. Lastly, in Sec.[V]we look 
forward to LBNE with its longer baseline as a way to break 
most of these degeneracies, and over most of the parameter 
space, deduce the hierarchy and presence or absence of NSI. 
Our conclusions are summarized in Sec.lVIl 



II. NSI IN THE e T SECTOR 
A. Standard Three-Flavor Vacuum Oscillation Pliysics 

Measurements from solar, atmospheric, reactor and beam 
neutrino experiments indicate flavor transformations between 
the three known flavors of neutrinos. Thus neutrinos are mas- 
sive and their mass and flavor bases are misaligned. In the 
mass basis, the three-flavor oscillation Hamiltonian is param- 
eterized as follows. 



^vac - 



-Aa-A 
Aq-A 
A0+Ay 



(1) 



where A = Amgj/ (4£'v), and Aq = Am2j/(4£'v), with m?- = 
mj — rrPj and the neutrino energy. We stress that the restric- 
tion to the known three flavors is a simplifying assumption, in 
keeping with our simplified framework. In general, one may 
want to consider NSI in the presence of additional sterile neu- 
trino states. 

In writing Eq. ([T]i we have used the fact that ^ m,, and 
dropped both the Majorana phases and an overall constant, 
since neither contribute to oscillation probabilities, as is well 
known. Although the magnitude of the atmospheric splitting 
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IAot^jI is well measured by MINOS EH to be (2.32+o og) x 
lO^^eV^, and the solar splitting Am^j = 7.1+^-^ x 10"^ eV^ 
is well measured by KamLAND 132), the sign of the former, 
which determines the ordering of the three mass states, re- 
mains unknown. We follow the prevailing convention in re- 



ferring to A > as the normal hierarchy (NH), and A < as 
the inverted hierarchy (IH). 

Eq.([I} can be rotated to the flavor basis i/™." UHIZ'V' 
via the PMNS matrix. Following the standard notation Ii33l . 



/I \ / C13 si3e-'^\ / C12 ^12 0\ 

UpMNS = C23 ^23 1 Sn Ci2 , 

\0 ~S23 C2ij V-^ne'^ C13 / V 1/ 



where Cij = cos0,-^ and Sij = sin0,y. Thus in vacuum neu- 
trino oscillations are characterized by three angles and one 
CP-violating phase 5. We adopt the values 023 = 45° ± 
7.7° ilB, 012 = 34° ± 1.1° 134|, and the recently updated 
013 = 8.7° ± 0.5° 135] (see ©I for the pubHshed results from 
Daya Bay). The CP-violating phase 5 remains unknown and, 
together with the type of mass hierarchy, constitutes the main 
"known" target for the next generation of the long-baseline 
experiments. As we will see below, NSI provide another im- 
portant target. 

B. Non-standard Oscillation Physics 

When neutrinos propagate in matter, their coherent forward 
scattering on the medium (refraction) must be taken into ac- 
count 1 1 1. The presence of NSI between neutrinos and quarks 
(and/or electrons) in the medium modifies this matter poten- 
tial, thereby altering the neutrino oscillation probability in 
long-baseline experiments. Consistent with our framework, 
we turn on only non-standard e — T interactions, 

/ 1 |e,.|e-'^v\ 
HZ^ V2Gf , (2) 

\\ee,\e'^^ 0/ 

where Gp is the Fermi constant, and n^, is the number den- 
sity of electrons. For the experiments considered here, the 
neutrino beams stay inside the continental crust, for which 
we take = lA mol/cm^. The "11" entry of Eq. (j2| is the 
SM piece, while gg^ parameterizes the strength of NSI in the 
electron-tau sector As akeady mentioned, the NSI piece can 



come from electrons, up or down quarks, or a combination 
thereof. We therefore define Set = 'T.f=e.u.d^ex{nf) /ng, such 
that the NSI couplings in Eq. (j2]i are normalized per elec- 
tron. As discussed in the introduction, the NSI Lagrangian 
term we have in mind is a effective operator of the form, 
-2\/2 e/i^Gf 7'' Vt) + /i.e. This is a useful pa- 

rameterization for our purposes since many possible UV com- 
pletions lead to such an operator at low energies. We leave the 
model-building and non-oscillation phenomenology of such 
UV completions to a future study. 

The rationale for choosing to consider only £ex is as fol- 
lows: (i) The interactions involving the muon neutrino are 
more strongly constrained than those of either Vg or v^. (ii) 
Small flavor-diagonal (FD) NSI for and/or Vf would lead to 
a small modification of the diagonal matter potential, as if the 
neutrinos travel through matter of slightly different densities, 
while small flavor-changing (EC) NSI can lead to qualitatively 
new effects, as will be seen later (iii) Large FD NSI could be 
allowed in combination with suitably chosen EC NSI lilOlfTTl . 
but as mentioned above, in this paper, for simplicity we do not 
consider cancellation between several NSI parameters. 

Observe that the NSI piece in Eq. (j2]) is in general complex, 
with a new irremovable CP violating phase 5v = arg [Eez) l3 
|37][38l. We shall see that that this phase interferes with the 
8 phase of the "atmospheric" oscillation term, producing im- 
portant effects and fundamental degeneracies. 

The measured value of 0i3 and the assumed smallness Eex 
leads us to anticipate the concomitant smallness of the prob- 
ability P (v^ Ve) which we will verify a posteori. With 
the Vg appearance probability assumed small we can treat the 
problem in perturbation theory, resulting in lfT2l : 



P v„ 



G\ sin 023 G2 cos 023 



Ai A2 
Gi = \/2GFn,|egT|e'^''cos023+Asin20i3e'^ 
G2 = V2Gf«,|e,^|e'^''sin023-A0sin20i2, 



(3) 

(4) 
(5) 



with A1.2 ~ A± |A| — V^Gfiie (neglecting the solar splitting A©, see below), and L the distance of propagation. Compared 
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FIG. 2: Here we examine MINOS sensitivity to NSI in the nor- 
mal (upper panel) and inverted hierarciiy (lower). The lighter re- 
gion comes from varying the vacuum phase 8 with SM physics only, 
while the darker regions come from including NSI with le^Tl = 0.4 
and varying the NSI phase with the vacuum phase set to zero. 

to lfT2ll we have taken only Igg^j nonzero and restored explic- 
itly the phases of both the vacuum and the NSI pieces. For 
typical energies Ey — 2 GeV, 623 = 7r/4, and 613 — 8.7°, the 
relevant parameters in the problem are 




Asin20i3 = 0.87 x 10"" eV, 

\/2Gf;vcos023 = 0.76xlO""eV, 
Ar^sin20,7 = 0.09 X 10"" eV. 



(6) 
(7) 
(8) 



The physics behind the general form of Eq. ^ can be 
understood as follows. The — ^ Ve conversion amplitude 
receives contributions from two frequencies, related to the 
"atmospheric" (Ai ~ 2A) and the smaller "matter" (A2 ~ 
— V^Cfne) splittings (the "solar" A© is smaller still and is 
for simplicity neglected). In the standard case (no NSI), the 
term Asin20i3e'^ drives the conversion Vfi Vg with the at- 
mospheric oscillation frequency, as captured by the Gi term 
in Eq. ([3]l. The smaller off-diagonal "solar" term A0sin20i2, 
captured by the G2 term, also drives the transition, but with a 
smaller frequency, A2. 

The standard CP violation search is based on the interfer- 
ence of the terms in Eqs. (j6]l and ([8]l. The magnitude of in- 



terference is dictated by the phase 5, which is responsible for 
CP violation, and by the oscillation phases, arg (e'^i 2i_ 
in the two channels. Furthermore, since the solar term ([8]l is 
an order of magnitude smaller than the atmospheric one, CP 
violation appears as a subleading effect, modifying the lead- 
ing probability due to Eq. ^ by at most ^ 20% (when the 
interference is maximally constructive or destructive). Its ob- 
servation thus requires sufficient experimental precision. 

The presence of nonzero Ect NSI modifies the amplitudes 
of both channels. Physically, ordinary oscillations generate 
Vt, which EeT then converts into Ve- With \eeT\ of order 0.2- 
0.4 and Ey ^ 2 GeV, the hierarchy of terms in Eqs. ( 6|7|8 1 
becomes: atm > NSI > sol. Thus, the expected NSI effect is 
still subleading, but in general larger than the standard signal 
of CP violation. The observable effect of NSI then depends at 
leading order on the relative phase 5y — 8. As an illustration, 
when this relative phase is zero and leerl = 0.2, one expects a 
^ 30% enhancement on top of the leading atmospheric prob- 
ability. 

We finish this section with two important corollaries to the 
above discussion. First, since ordinary oscillations form the 
necessary first stage of the conversion, at high neutrino energy 
Ey, the Ve —> Vfj conversion probability goes to zero even in 
the presence of nonzero Eer- Thus, while naively one might 
expect nonstandard matter effects to be cleanly manifested 
at high energies', this is not so for Eer- The best energy to 
probe EeT in Ve conversion is at the appearance maxi- 

mum. Thus, NOvA (and and its proposed successor, LBNE) 
are suitable experiments to look for this type of new physics. 

The second observation is that the recent measurement of 
large On is crucial in giving MINOS, NOvA, and LBNE sen- 
sitivity to NSI, since the NSI-driven conversion interferes with 
the "standard" amplitude driven by 613. As an illustration, 
consider the fact that with 613 = and leg^l = 0.2 this prob- 
ability at NOvA is below 0.005, even with constructive NSI- 
solar interference. This signal is certainly below the sensitiv- 
ity reach of the next generation of long-baseline experiments. 
To get an observable signal at MINOS, P (v^ V^) 0.05, 
with (?i3 = requires large NSI, {Eerl ~ 0.9 [13 ■ Since in 
the past year the value of 613 was measured to be sufficiently 
large, it is time to revisit the sensitivity of MINOS to NSI. 



III. MINOS 

Of the long-baseline oscillation experiments that already 
have data, MINOS provides the best sensitivity to NSI. This is 
due to their relatively long baseline and the resolution of their 
P{Vii Ve) measurements. In fact, as we show below, the 
Ve appearance search by MINOS [39J has already started ap- 
proaching the region of the parameter space favored by solar 
data. 

We begin by asking what the NSI sensitivity of MINOS is in 



Which is indeed true for, e.g., 
conversion at high energy. 



In that case, one expects 
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the case of the normal hierarchy, A > 0. To this end, Eq. ^ is 
plotted for the MINOS baseline in the top panel of Fig.|2] both 
with NSI of magnitude je^.^ | = 0.4 (darker region) and without 
it (lighter region). The standard physics region is swept out by 
the variation of the vacuum phase 5, while the NSI region has 
vanishing vacuum phase but varying the 5v phase. All other 
oscillation parameters are set to their measured central values 
(see Sec. II A i. From Fig. |2] we see that in the presence of the 
NSI the two hierarchies overlap considerably. This is the first 
example, of many to be encountered, illustrating the ability of 
NSI to easily confuse the interpretation of experimental data 
when viewed in terms of standard oscillation physics. 

Furthermore notice that when NSI interferes constructively 
with the SM contribution, the — > Ve conversion probabil- 
ity can be as large as ^ 0.11. At face value, this is in conflict 
with the 90% CL bound from MINOS |39|. However, the sen- 
sitivity of MINOS diminishes appreciably after marginalizing 
over the uncertainty in known oscillation parameters. Indeed, 
recall that Fig.[2]does not include the allowed variation in the 
standard oscillation parameters, 613 and 023- With the now 
small error bars on 613 from Daya Bay, tuning 613 down does 
not change the conversion probability appreciably. Suppose, 
however, that the actual 023 value in nature is lc7 away from 
its central value, 023 ~ 37°. This renders the largest NSI prob- 
ability P (v^ Ve) « 0.083. This conversion probability is 
only 1 (7 higher than the MINOS best-fit point in the NH. A 
similar exercise in the case of the IH does not yield any exclu- 
sion of the solar-preferred values of NSI. 

Notice, additionally, that the Vg appearance probability in a 
given hierarchy can mimic that of the other hierarchy. Com- 
paring the two panels of Fig.|2]one observes that the presence 
of NSI significantly exacerbates the difficulty in determining 
the hierarchy. Folding in information regarding the location 
of the peak would certainly increase the sensitivity of MINOS 
to the sign of the hierarchy, since the NH peaks at somewhat 
lower energies. Given their substantial background however, 
MINOS appears unlikely to have sufficient sensitivity to this 
signature of the hierarchy. 

A priori, MINOS could have observed indications of 
anomalous conversion due to NSI \eeT: \ ^ 0.4, had the phases 
been such that the interference is constructive and the true 
value of 023 were on the high end of the allowed range. This 
shows that MINOS is approaching an interesting range of sen- 
sitivities. With MINOS already having skimmed off a small 
portion of the NSI parameter space thanks to large 0i3, it is 
natural to ask how this progress will be improved upon. The 
next relevant experiment is NOvA, which is designed to im- 
prove upon the MINOS — > Ve results by about an order of 
magnitude. Thus, it is timely to examine the effect of the NSI 
at NOvA. 



IV. NOvA 

With a baseline of 810 km, a narrow-band 2 GeV beam 
around the expected v^ — >■ Vg conversion maximum and en- 
hanced detector sensitivity, the NOvA experiment |40| will 
offer the next opportunity to measure or constrain NSI. Sup- 
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FIG. 3: Here we have fixed the neutrino energy to 2 GeV, and plot- 
ted the ensuing values of /'(v^ Ve) and P{Vfi — )■ Ve) for NOvA. 
The solid ellipses refer to SM only interactions, for the normal (blue) 
and inverted (red) hierarchy, with the vacuum phase 8 varying con- 
tinuously along each of the ellipses. For the dashed ellipses we have 
used the NSI values \£er\ = 0.4 and 8v = 0, corresponding to the 
solar neutrino fit in Fig.[T](solid curve there). 



posing NOvA runs for 3 years in each mode, we ask what 
additional information can be gleaned from combining data 
from both neutrinos in P{Vn — > Ve) = P and antineutrinos in 
P{v^, ^ V,) = p. 

A particularly simple way of depicting this is through the 
use of a "bi-probability plot" 141]. This method has been 
widely used by the NOvA collaboration for presenting their 
standard physics expectations. With only SM interactions, a 
bi-probability curve is constructed by fixing the neutrino en- 
ergy Ev and the mass hierarchy and varying the vacuum phase 
5. This "theory only" prediction can be used to understand the 
required experimental sensitivity to the phase 5 and the sign 
of A. Notice that although a bi-probability plot only includes 
information from a single energy bin, for NOvA this is not a 
significant limitation since it is by design a narrow-beam ex- 
periment with Eyr^l GeV (cf. Fig. 4.5 in f40l). 

The predicted bi-probability curves for the NOvA baseline 
are depicted in Fig. [3] The solid curves assume standard os- 
cillations physics only, with the blue ellipse referring to the 
NH and the red ellipse to the IH. The perimeter of each curve 
is traced from the variation of the vacuum phase, < 5 < In. 
The points closest to the center of each ellipse correspond to 
no CP violation (5 or n), while those farthest from the 
center correspond to maximal violation (5 — 7t/2 or 3;r/2), 
cf. Fig. 3.3 in 1 40 1 . Notice that as 5 runs through its range, the 
probabilities are modulated by at most ^ 0.01, or about 20%. 
This is indeed expected from the interference of the solar term 
with the leading atmospheric term, as discussed in Sect. |IIB| 

This "theory only" plot could in fact serve as a guide for the 
experimental sensitivity, if combined with the expected reso- 
lution information. The collaboration projects its 2(7 resolu- 
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tion to be roughly a circle in the bi-probability plane, with the 
diameter of ^ 0. 15 [42J. For example, by measuring P ^ 0.05 
and P ^ 0.03, NOvA could claim to have evidence that the hi- 
erarchy is normal at almost 3(7 confidence level (C.L.). More- 
over, at 2 (J C.L., the phase 8 would be determined to be 
3n/2±n/2. As a second example, if NOvA were to mea- 
sure P ^ P ^ 0.04, the hierarchy could not be established. It 
would be known, however, that the 5 phase is either in the in- 
terval 7r/2±7r/2with NH, or in the interval 3;r/2±;r/2 with 
NH. In other words, the existence of CP violation would be 
suggested at the 2a C.L. 

Crucially, all these determinations would apply only if it 
were somehow known that NSI was ruled out. When we allow 
for nonzero NSI in the e — T sector, the situation changes con- 
siderably, and the regions corresponding to the two hierarchies 
expand significantly. To illustrate this, we first consider the 
NSI scenario that was used in the fit of the solar data in Fig.[T] 
Recall that the thick curve in that figure was constructed for 
ICerl = 0.4 and the phase 5v = 0. In this case, at NOvA the 
variation of the vacuum phase traces out the dashed ellipses in 
Fig. [3] analogously to the SM ellipses discussed above. When 
NSTSM interference is constructive the probabilities in both 
neutrino and antineutrino modes can be substantially larger, 
approaching 0.08 — 0. 10. If nature prefers such fortuitous val- 
ues of NSI, NOvA will see dramatic deviations from the SM 
expectations. We see that NOvA is in an excellent position to 
probe NSI couplings at a level suggested by the solar data. 

Notice that the size of the NSI effect in Fig.[3]can be once 
again roughly understood from the discussion in Sect. |IIB| 
since ^/2Gfng cos 023£eT is ^ 3 times larger than A0 sm2di2, 
the dashed ellipses are roughly a factor of three longer 

Alternatively, one may be interested in the general question 
of NOvA's sensitivity to Eez, without reference to the solar 
data. To address this issue, we vary the 5v phase in its entire 
range. We also restrict \£eT\ to a smaller value, 0.2, which is 
justified given the level of sensitivity seen in Fig. [3] With this 
new fiducial value, we can repeat the same exercise as before 
but now with several different values of the v-phase. Doing 
so in Fig.|4] we find the regions shown in Fig. |4] (fop). 

Doing this for all possible v-phases from to 2n traces out 
the large blobs for the two hierarchies, see the bottom panel of 
Fig.|4] Note that for any given point in the NSI regions, many 
ellipses with distinct vacuum and v-phases intersect. This 
means considerable degeneracy of parameters corresponding 
to distinct values of both the v-phase and the vacuum phase 8 
intersecting at any given point. This is especially true in the 
central regions where one can see what in mathematical jargon 
is known as a double "swallowtail catastrophe," as the param- 
eter region folds on itself. Based on this observation alone, 
NOvA will have difficulty in knowing what combination of 
phases they are truly measuring. 

Let us now consider four qualitative possibilities for the 
outcome of NOvA as illustrated in Fig.|4](/7offom): 

(1) Clear NSI and hierarchy determination. For example, 
the point = (0.08,0.04) (o in Fig. [iji is widely dis- 
crepant from any SM-only explanation. Moreover, such 
a set of probabilities would not be likely to have come 
from SM or NSI cases in the inverted hierarchy. In this 
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FIG. 4: Similar to Fig. [5] although here we fix \Eet\ = 0.2. Top: 
Scanning over many values of the NSI phase Sy results in signif- 
icant degeneracy, with many distinct curves intersecting at a given 
point. Bottom: The result of allowing 8v to float arbitrarily yields 
the shaded blobs. See text for a discussion of the four qualitatively 
distinct points. 



case, one could fairly say that the normal hierarchy and 
nonzero NSI are strongly favored. This represents an 
example of the best case scenario for NOvA being sen- 
sitive to this new physics. 

(2) Clear NSI determination only. Here the point {P,P) = 
(0.055,0.06) (• in Fig.|4|, would indicate a strong pref- 
erence for the existence of NSI. Of course with this 
measurement alone however, no confident statement 
about the sign of the mass hierarchy could be made. 

(3) Hierachy determination only. An exemplary point of 
this possibihty is offered by (f,P) = (0.06,0.03) (■ in 
Fig.|4]). Here the normal hierarchy would be mildly pre- 
ferred over any explanation based on the inverted hier- 
archy. However, this point is degenerate with SM and 



7 



Degeneracies at NOvA, v mode 



0.00 



t 
1 




S=7T \e^J=0 


1 




(5=3.98, |e„|=0 .2. i5,,=;r/4 


1 

_ 1 




5=0.55, |e„l=0.2, 5,-11 


1 1 
11 
11 


/ /'""^^^ 


5=2.14, |eerl=0.2,5v=3;r/4 


1 1 
1 1 
1 1 

1> ' 
1' / 
\\ 11 


* / ... 

w ' 

'/ ' •*lk. 
1 ' .• 

1 •' ^ 

1 / 
/ • 




. i 

\\\ ! 
U' .7 / 


/ / 
/ / 

• 




K .7/ 
■ft ' ' 













0.5 1.0 1.5 2.0 2.5 3.0 3.5 4.0 

Ey (GeV) 

Degeneracies at NOvA, v mode 



t 0.04 
0.02 
0.00 





6=71, \e„\=0 


i'- 

•1 '. 


5=3.98, |e„l=0.2, 5„=ff/4 


- ; 1 '; 
■ 1 ; 
1 ■ 

1 ! 

\ 1 '. 
I 1 ■< 


5=0.55, |eerN0.2,5v=l.l 
5=2.14, |eeTl=0.2, 5„=3;r/4 


I '* ' 




'.li /// / 




111*. // / / 

'i\ •■'7' .' 





0.5 1.0 1.5 2.0 2.5 3.0 3.5 4.0 



(GeV) 

FIG. 5: An illustration of some of the NOvA degeneracies at 
2 GeV in v mode {top) and v mode {bottom). All curves as- 
sume the normal hierarchy. The blue (solid) curve is without NSI 
and setting the vacuum phase to 5 = ;r. The remaining curves all 
have NSI of the same magnitude \£i,x\ = 0.2, but take (S,8y) = 
(3.98, n/A), (0.55, 1.1), (2.14,3;r/4), in the green (dashed), red (dot- 
ted), and magenta (dash-dotted) curves respectively. Note that the 
choice of phases produce degenerate results at 2 GeV in both neutri- 
nos and antineutrinos. 

NSI explanations. 

(4) No NSI or hierarchy determination. The point — 
(0.04,0.03) (▲ in Fig|4]) is an example of one of the 
worst cases for NOvA to have a clear signal of any as 
of yet unknown parameters. At such a point, one cannot 
rule out the existence of NSI or establish the sign of the 
hierarchy. 

We further illustrate the degeneracy of the last point ▲ in 
Fig.|5] where we plot the conversion probability as a function 
of energy. One of the curves has only standard physics and 
5 — n, while the other three curves have NSI with different 
choices of the vacuum and v-phase. Both neutrino and an- 
tineutrino curves intersect ?A Ey —2 GeV and with a naiTow 
band setup of NOvA give the same event rates. Importantly, 
the degeneracy is not absolute and can be broken by measure- 
ments at different energies and/or baselines. We will return to 
this in the next section. 



NOvA, Vary 6*23 
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FIG. 6: Here we have fixed the neutrino energy to 2 GeV, and plot- 
ted the ensuing values of /'(v^ ) and P{v^ ) for NOvA. 
The outer (inner) cylinder regions refer to SM only interactions, for 
the normal (blue) and inverted (red) hierarchy, the angle 623 vary- 
ing within its presently (large, light cylinders) and future (small, 
dark cylinders) allowed 90% CL region centered on n/A. The larger 
shaded regions come from fixing \eeT\ = 0.2 and varying both the 
vacuum and the matter phases. 



Before concluding our discussion of NOvA, it is worth 
mentioning yet another type of degeneracy, which exists be- 
tween the standard oscillation parameter 623 and NSI. The un- 
derlying physics behind this degeneracy is evident from the 
analytical form of the probability in Eq. ([3]i: the change of G\ 
due to NSI can be partially undone^ by appropriately modify- 
ing the factor of sin 2 023- We have akeady seen an example 
of this degeneracy in the case of MINOS, where the effect of 
NSI could be partially undone by adjusting the value of 023 • 
Let us now describe, quantitatively, this degeneracy at NOvA. 

Once again, we turn to the bi-probability plane. The ef- 
fect of varying 023 in this plane - assuming standard physics 
only - is to shift the solid ellipses in Fig.|3]toward and away 
from the origin {cf. |42|). In other words, varying 023 turns 
the ellipses into cylinders. We illustrate this in Fig. |6] which 
generalizes Fig.|4]to the case of uncertain 023. The filled back- 
ground regions are once again obtained by varying 5 and 5v, 
assuming fixed \£ex \ — 0.2 and 023 = n/A. The lightly shaded 
foreground cylinders are the result of setting \eeT\ to zero and 
varying 8 over its full range and 023 over the range allowed 
by the Super-Kamiokande atmospheric data^ sin^ 2023 > 0.93. 



^ The degeneracy is clearly partial, even for unconstrained 623, since sin 2 623 

is real, while G[ is complex. 
^ The values of sin^ 2023 infeired from the Super-Kamiokande data actually 
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We see that the degeneracy between NSI and 023, while only 
partial, reduce the possible statistical significance of the NSI 
signal. 

It then becomes crucial for the NSI search to have a better 
determination of 023- Fortunately, NOvA itself will be able 
to improve the error on 023 by about a factor of two, by us- 
ing data in the P(v^ V^) channel |42|. This results in the 
smaller darker cylinders also seen in Fig. |6] which correspond 
to sin'^2023 > 0.98. Thus, the improved measurement of 023 
by NOvA translates into improved NSI sensivity. 

Further improvements could be achieved in two ways: by 
tightening the bounds on 023 from other measurements and by 
going to a longer baseline. The latter possibility is described 
in the next section. As for the former avenue, an important 
contribution could be made by the Deep Core component of 
IceCube. The authors of |43| estimate that with 10 years of 
data. Deep Core will be able to measure 023 to within ± 2" 
from maximal mixing. This would be an enormous help in 
conjunction with the data from NOvA in measuring or con- 
straining NSI and the sign of the hierarchy. 

The flip side of the 023— NSI degeneracy is that under the 
assumption of NSI, NOvA loses any power of "octant" dis- 
crimination (i.e., discrimination between the light side, 023 < 
;r/4, and the dark side, 023 > ;r/4, in the terminology of |44|). 
For example, a point like {P,P) = (0.06,0.04), which under 
the standard assumptions could be taken as an indication that 
023 is in the dark side, could instead be coming from NSI and 
023 on the light side. 

Of course, uncertainties on the precise value of 0i3 remain 
as well, though these do not produce effects nearly as big as 
the uncertainties on 023. This is due to the fact that 0i3 is 
known to within 6% lf35l already, whereas 023 is known only 
within 17%. Importantly, since the Daya Bay baseline is short, 
its measurement of 0i3 is not impacted by the existence of 
NSI. 



V. BREAKING NOvA DEGENERACIES WITH LBNE 

The Long-Baseline Neutrino Experiment (LBNE) l! 45l will 
be a very useful tool in breaking many of the degeneracies 
that could easily exist after NOvA. In this section we illus- 
trate how this can be done with the Fermilab-to-Homestake 
baseline (1300 km), as it yields the most sensitive probe of 
matter effects of the proposed LBNE baselines. 

Before considering NSI, let us first examine LBNE's SM- 
only sensitivity to the hierarchy. While LBNE is not planned 
to run with a narrow-band beam, to facilitate comparison with 
our NOvA results we begin by showing the predictions in the 
bi-probability plane for Ey —2 GeV. The result is depicted in 
Fig.|7]and is to be compared with Fig.|6] 
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change for large NSI, as shown in IIOIIIII We neglect this change here, 
since we work in the assumption of small NSI (no cancellations between 
different e's). 



FIG. 7: Top: Here we have fixed the neutrino energy to 2 GeV, and 
plotted the ensuing values of /'(v^ — > Vf ) and -P(v^ — > v?) for the 
Fermilab-to-Homestake baseline, 1300 km. The foreground cylin- 
ders refer to SM only interactions, for the normal (blue) and inverted 
(red) hierarchy, with the vacuum phase and 623 varied as in Fig. |6] 
The larger shaded regions come from fixing \£i,^\ = 0.2 and varying 
both the vacuum and matter phases. 



As can be seen from Fig. |7] the cylinders for the two hi- 
erarchies (obtained by varying 023 as in the last section) are 
very well separated. Thus, under the assumption of vanishing 
NSI, the Homestake baseline has the capability to definitively 
determine the sign of the hierarchy. With the addition of the 
NSI, the overlap between the two regions is restored. Thus, 
for part of the parameter space, the fundamental ambiguity in 
determining the hierarchy remains. Yet, a measurement of, 
say, {P,P) = (0.06,0.06), would actually be a very exciting 
development, as it would indicate the existence of NSI. 

Notice that for a significant part of the parameter space, 
the presence of NSI would be even more apparent. For ex- 
ample, notice in Fig. [7] that the most extreme NSI point in 
the normal hierarchy is more than 0.05 away from any SM 
physics interpretation. Thus if LBNE can achieve a sensitiv- 
ity of f ^ 0.01 to Ve appearance they could observe a near 5g 
detection of nonzero NSI. Furthermore, at this extreme point 
one would have some information on the magnitude of NSI 
since the P(v^ Ve) is never this large with Set =0.1. In- 
deed LBNE is so sensitive to NSI that it may be able to have 
a 3(7 signal of NSI even with Set = 0.1. 

Let us now consider the full spectral information and reex- 
amine the example of degeneracy at NOvA depicted in Fig. [5] 
Plotting the same set of curves as before, but now at the longer 
Homestake baseline we see, in Fig. [8] that the four distinct 
possibilities become much easier to discriminate, even with 
only information at 2 GeV. Beginning in the v mode, we see 
that some mild ambiguity continues to exist in discerning be- 
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FIG. 8: An illustration how LBNE can break the NOvA degenera- 
cies depicted in the bottom panel of Fig. |4] Here the hierarchy is 
normal and the color scheme is the same as in Fig.|4] 



tween the red curve with {5,5y) — (3.98, 7t/4) and the green 
curve with {5,5y) = (0.55, 1.1). This is an excellent exam- 
ple, however, of the benefit of having broader energy sensitiv- 
ity. If one has access to low-energy information, especially as 
low as the second oscillation maximum, then one could eas- 
ily tease out the difference between the two curves. Likewise, 
information in the antineutrino mode (bottom panel of Fig. |5]l 
demonstrates an orthogonal method for the discrimination of 
the phase. Notice that in the v mode, the two closest curves in 
the V mode are the farthest apart. 

We see that the data from NOvA and LBNE should be an- 
alyzed together. In fact, the combined analysis of data from 
several baselines, such as T2K, NOvA, LBNE, and Deep Core 
promises to provide a very powerful search tool for NSI. 



VI. CONCLUSIONS 

We have argued that present solar neutrino data may pro- 
vide a hint of new interactions between matter and neutri- 



nos. Such interactions are not well-bounded by existing ex- 
perimental probes. In this paper, we have advocated a delib- 
erately simplified, effective framework for gauging the sensi- 
tivity of long-baseline neutrino experiments to NSI. We have 
further argued that the recent measurements from Daya Bay of 
large 613 magnify the effects of NSI via quantum interference 
with standard vacuum oscillation physics. We have found that 
near-term long-baseline experiments will be an important test- 
ing ground for the viability of NSI. In particular, the v^ — ?► Ve 
search by MINOS is already starting to exclude previously al- 
lowed regions of parameter space. The upcoming NOvA ex- 
periment will have increased v^, — > Ve sensitivity, being able 
to see 3(7 deviations from the no-NSI hypothesis. However, 
the additional Cf -violating phase 5v from NSI complicates 
the determination of the neutrino mass hierarchy, the octant 
of 023 and the vacuum phase 5. We discuss how this inher- 
ent degeneracy could be broken with information from other 
baselines. To this end, we demonstrated that a long-baseline 
experiment with L ^ 1300 km (LBNE) would be able to break 
many of these degeneracies and determine the hierarchy and 
the presence of NSI. Combining data from other baselines is 
expected to resolve the degeneracies even further. 

Lastly, we note while we only considered the impact of 
new physics on neutrino propagation in matter, depending on 
the underlying physics model of NSI, one may also encounter 
other new physics effects. In particular, if a new light sector 
is responsible for the discrepancy in the solar data, it could 
also play an important role in neutrino production and detec- 
tion (see for example L271), though existing data from fixed- 
target experiments Il46l and so-called "zero distance" oscil- 
lation experiments ^47\ provide strong constraints on some 
models of such new physics. We note that distinct physical ef- 
fects on neutrino production, detection, and propagation could 
be discriminated with a long-baseline experiment if far detec- 
tor data is used in conjunction with data from a near detector. 
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